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^Nj ■ Abstract. We consider a Laplace operator on a random graph consisting of infinitely 

many loops joined symmetrically by intervals of unit length. The arc lengths of the loops 
(j^ are considered to be independent, identically distributed random variables. The integrated 

■ density of states of this Laplace operator is shown to have discontinuities provided that the 

distribution of arc lengths of the loops has a nontrivial pure point part. Some numerical 
illustrations are also presented. 
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^ ! 1. Introduction 

m . 

O ■ In this paper we study a model which perhaps provides the simplest example of a differ- 

ential operator on a nontrivial infinite random metric graph. The graph consists of infinitely 
^T) I many loops joined symmetrically by intervals of unit length. The arc lengths of the loops 

O ■ are independent, identically distributed random variables. A similar model where the arc 

lengths of the loops are kept fixed, was considered by Avron, Exner, and Last in [2]. They 
Oh' called this model a Necklace of Rings. Mimicking this terminology we will use the name 

. Random Necklace for the model considered here. 

I The main objective of the present work is to study the integrated density of states and the 

Lyapunov exponent for the Random Necklace Model. It is well known that the integrated 
density of states for metrically transitive Schrodinger operators on the line is continuous 
at all energies since the multiplicity of their spectrum is not greater than two. In contrast, 
^ ■ the Laplacian of the Random Necklace Model can have eigenvalues of infinite multiplicity. 

Therefore, the integrated density of states may have discontinuities at those energies. We 
will show that this is indeed the case provided that the distribution of arc lengths of the loops 
has a nontrivial pure point part. An explicit description of the set of all energies where the 
integrated density of states is discontinuous will also be given. Moreover, we will show that 
the perturbation of the Laplacian by a magnetic field in general smoothes out the integrated 
density of states such that some of its discontinuities disappear. 

Discontinuities of the integrated density of states for some discrete random Laplacians on 
Delone sets have been studied recently by Klassert, Lenz, and StoUmann in [9]. Earlier such 
discontinuities had been observed in [1], [7], [10], [18] for discrete Laplacians associated 
with Penrose tilings. The appearance of discontinuities is again related to the existence of 
infinitely degenerated eigenvalues. 

The plan of the present work is as follows. The model is defined in Section 2. In Section 
3 we decompose the integrated density of states N{E) into the integrated density of loop 
states A^'°°P(i?) (i.e., the states supported on loops of the graph) and the remainder N{E). 
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Next we prove that A^'°°p [E) has discontinuities if and only if the distribution measure has a 
nontrivial pure point part. In Section 4 we show the smoothing effect of a constant magnetic 
field. In Sections 5 and 6 we accommodate the scattering theoretic method proposed in [13] 
(see also [12]) to calculate the integrated density of states and the Lyapunov exponent for 
the Random Necklace Model. Using this approach we prove the positivity of the Lyapunov 
exponent for almost all energies and study the set of its zeroes. In Section 7 we prove the 
continuity of N{E) and discuss its further regularity properties. 

We present also some numerical illustrations performed by the method developed in [13]. 
Part of the material presented here has previously appeared in [1 1]. 

2. Random Necklace 

In this section we give a precise formulation of the Random Necklace Model and discuss 
some of its basic properties. Consider an infinite graph F consisting of loops Lj with arc 
lengths 2uj joined symmetrically by unit intervals Ij = [0, 1] (see Fig. 1). Any loop will 

be realized as a union of its upper L^^^ = [0, uj] and lower L^~^ = [0, uj] parts. We 
always assume that the left vertex of any loop Lj corresponds to the point x = for all 

three bonds Lj^\ Lj \ and Ij adjacent to this vertex. Thus, the right vertex of the loop Lj 
corresponds to the point x = 1 of the interval /j+i and to the point x = uoj of the intervals 

L^^\ Further, we suppose that lo = {LOj}j^z forms an i.i.d. sequence of random variables 
with distribution measure x and satisfying < Uj < K. The underlying probability space 

is, therefore, Vl = [0, K]^ with the product probability measure P = Xj-g^x. 
With the graph T we associate the Hilbert space %, 

(I) 'K = ^L\lj)®L\Lj) with L\Lj)=L\Lf^)®L\L^r)). 
According to the decomposition (1) we will write the elements V' of IK as follows 

where i^f G L'^{Ij) and ^ e L'^{Lf^)- 



illustration.png 



Fig. 1. Random necklace. 

Note that the Hilbert space "K depends on the random variable uj = {uoj}j^i. By proper 
scaling it is possible to formulate the problem equivalently on a Hilbert space independent 
of UJ and work with a random operator on a deterministic graph. However, we will not use 
this formulation. 
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On 3i we consider the negative Laplacian — A(ti;) as the operator of second derivative 
with local boundary conditions of the form 



/V'f (0) 



4+)(o) 

(-), 



(2) A 
at the left vertex of the loop Lj and 



+ B 



0, 



(3) 



A 



+ B 



( ~^f{D \ 



0, 



at the right vertex. Here A and B are such that they define so called standard boundary 
conditions, 



-1 



(4) 



A 



1 




and 



B 








By a general result in [14] the operator — A(ti;) is self-adjoint and nonnegative. By the 
general theory of metrically transitive operators (see, e.g., [22]) the spectrum of — A(a;) as 
well as its components are deterministic. 



It is easy to see that the numbers Eji 
with the eigenfunctions 



(nl/uij) , j,l G Z are eigenvalues of — A(a)) 



(5) 



± sm and 



for all k G 



Each of these eigenfunctions is compactly supported on the j-th loop. Therefore, we will 
call these eigenvalues Eji the loop eigenvalues. In fact, — A(u;) has no other eigenvalues 
with compactly supported eigenfunctions. The fact that there are no other eigenvalues with 
eigenfunctions supported on a single loop can be verified directly. Moreover, in Section 5 
below we will prove the following proposition. 

Proposition 2.1. Let tp be an arbitrary solution of the Schrodinger equation = E^. 
Assume that ^ on an open subset of some interval Ij. Then in almost all 
internal points of every interval 1^, A; G Z. 



3. Discontinuities of the Integrated Density of States 
Let r™'" with m, n G Z, — m < n be the graph consisting of n + m + 1 loops Lj of arc 



length uoj, j 



-m. 



, n joined symmetrically by intervals of unit length. Let — A*"'" be 



minus the Laplacian for the graph on T"*'" with the boundary condition (4) at all vertices 
except those vertices of the loops L-m and L„, where no intervals are attached. At those 
vertices we impose Dirichlet boundary conditions. 

Obviously, the operator — A"^'" has a discrete spectrum. Therefore the finite volume 
integrated density of states is well-defined: 



tr E_A'">"((— oo, E)) 
m + n + 1 
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where E_a'"."((5) denotes the spectral projection of the operator —A™'" corresponding to a 
Borel subset (5 C M. By standard arguments (see, e.g., [8]) one can prove that the Umit 

lim A^"^'"(S) =:N{E) 



m,n— >oo 



exists almost surely for all points of continuity of N{E). This limit is deterministic and is 
called the integrated density of states. At possible (at most countably many) discontinuity 
points of N{E) we make the convention N{E) = N{E — 0). 

Let P denote the orthogonal projection in onto the subspace generated by the eigen- 
functions corresponding to loop eigenvalues. Then we can decompose the integrated density 
of states into two parts 

trPE_^,„.„^^,j((-^.i?)) 



(6) 
and 
(7) 



iVl°°P(E) = lim 

m,n—» oo 



N{E) 



lim 



m + n + 1 

tr (/-P)E_A^.n(^)((-oo,£;)) 
m + n + 1 



such that N{E) = N{E) + iVi°°P(£;). It is quite easy to calculate N^°°p{E) explicitly 
Indeed, for £^ > we have 



iV^°°P(E) = lim (m + n + 1)"^ V #{Eji\ Eji < E} 

UJj^/E 



j=-m 
n 



lim (m + n + 1) ^ > 



m,n—> OO 



TT 



where \t \ denotes the largest integer strictly smaller than t, \t\ < t. By the Birkhoff ergodic 
theorem the Umit exists almost surely and is equal to 



(8) 



N^°°P{E) = / 
Jo 



uooVe 



TT 



dx(a;o). 



From (8) it follows that N^''°p{E) = for all E < tt^/K^, where K is the suppremum of 
the topological support of the measure x. 

We have the following elementary result on the integrated density of loop states. 

Proposition 3.1. Let E >Tr'^/K^. For all sufficiently small e > 

N^°°P{E + e) - N^°°P{E) = x{Se), 

where 

Se= U 



k=l 



VETe' ^/E_ 



and lt\ denotes the largest integer not exceeding t, \t\ < t. 
The Lebesgue measure of the set is obviously bounded by 

[KVE/7t\ 

-^^<-il+KVEM. 

k=l 
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Therefore, Proposition 3. 1 implies that if the measure x is purely continuous, then A^'°°p(£') 
is continuous. Moreover, if x is purely absolutely continuous with bounded density, then 
jY'°°P (E) is Lipschitz continuous. On the other hand, one can construct purely absolutely 
continuous x's with unbounded density such that the integrated density of loop states is not 
Holder continuous of any prescribed order. For instance, fix some a G (0, 1). If (ix(wo) = 
^jCjlwo — ii|~°di^o with ti being dense in {0,K) and Cj > satisfying Cj < oo, 
then N^°"^{E) is Holder continuous of order 1 — a and not of order /3 > 1 — a. Also, 
there are singular continuous x's such that the integrated density of loop states is not Holder 
continuous of any order a > 0. 

Proof of Proposition 3.1. For given E > choose e > so small that 



(9) 
Then 



and therefore 
and only if 



K 

TT 









TT 




TT 



( V^B+i - \/e) < 1 
< 1 for any < uq < K. The equality sign occurs if 



UJQ^/E + £ LUo\^ 

> k> 



TT 



TT 



for some integer k. Thus, 



'ojQyjE + e" 




ujQ\fE 


TT 




TT 



X5e('^0), 



where xSe is the characteristic function of the set S^. The claim follows from equation 
(8). □ 

Theorem 3.2. Assume that the probability measure x has a nontrivial pure point part Xpp, 

oo oo 

(10) xpp(-) =^Pi(5s,(-) with pi>0, 0<^pi<l, 

j=l i=l 

where Ss^ denotes the Dirac measure concentrated at Si. Consider the nonempty set 

\2 



(11) 



D^:={E eR+\E = {Trk/siy, Siy^O, piy^O, keN}. 



Then is the set of all points of discontinuity of the integrated density of states N{E) and 
for any E G 



(12) 



N{E + 0) - N{E) = am > 



i=l 



with ai = 1 if SiV E ji: is integer and = otherwise. 
Proof. From Proposition 3.1 it follows that 

iVi°°P(E + 0) - Ar^°°P(£;) = x(M), 
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where M = {ujq\ lvoVE/tt is integer}. Obviously, for any E the set M is at most discrete. 
Therefore, we have 

oo 

(13) iVi°°P(E + 0) - N'°°P{E) = Xpp(M) = ^ aiPi ^ 

i=l 

if E & and zero otherwise. 

Since N{E) is non-decreasing it follows that N{E) is discontinuous on D^. Moreover, 
we obtain (12) with "=" being replaced by ">". Thus, to complete the proof it suffices to 
show that N{E) defined by (7) is continuous. We will prove this fact in Section 7. □ 

Since the set M is discrete the sum in equation (12) involves only a finite number of 
nonvanishing terms. 

If the distribution measure x has a nontrivial singular continuous component, then the 
canonical decomposition of N^°"^{E) contains a nontrivial singular continuous part. This 
is established in the following proposition. 

Proposition 3.3. Ifx is purely singular continuous, then N^°°p{E) is a singular continuous 
junction. 

Proof. Since N^°°P{E) is non-decreasing, it is differentiable for Lebesgue almost all G M. 
By Theorem 3.2 it is continuous. Therefore, to prove the claim it suffices to show that 

dN^°°P{E)/dE = for a.e. E. 

Fix an arbitrary Eq > jK'^ and consider 

/cvr fcvr 



Fkie) = X 



keN 



as a function of £ G (0, r(£'o)) with r{Eo) being defined in (9). Since x is singular contin- 
uous, we have Fl^{e) = for a.e. e. From Proposition 3.1 it follows that 

7Vi°°P {Eo + £) = 7V'°°P {Eo)+ J2 

k=l 

for all £ G (0, r{Eo)). Thus, for almost all E G {Eq, Eq + r{Eo)) we have 

dN'--p{E) dAr'°op(£;o + £) L^^/^J 

dE = de = Fk{^) = 0- 

k=l 

This proves the claim since Eq is arbitrary. □ 

4. Magnetic Field Smoothes the Integrated Density of States 

Without loss of generahty we can assume that the graph T is imbedded in a plane in 
and the loops Lj are circles. In this section we consider a magnetic perturbation — A(a;; S) 
of the Laplacian described in Section 2. Assume there is a constant magnetic field S per- 
pendicular to the plane containing the graph. Since the j-th loop encloses an area i^I/tt the 
magnetic flux through this loop is given by = uj'j'B/'K. As shown in [17] prescribing 
magnetic fluxes through all loops of the graph defines a magnetic Laplacian uniquely up to 
a gauge transformation. The resulting magnetic Laplacian is again defined to be the operator 
of second derivative but now with different boundary conditions. The boundary condition 
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(2) at the left vertex of the loop Lj remains unchanged and the boundary condition (3) at the 
right vertex takes the form 



A, 




I 



with 
(14) 







A, 



'\ -1 

e'*^/2 -e-^*j/2 



,0 







-^f '(1) ^ 













1 g-j*,/2 



It is easy to see that the numbers Eji = {tvI/uj)^ are eigenvalues of — A(a;; "B) if and 
only if 'Bujj /ir"^ is integer. The corresponding eigenf unctions are given by (5). Also, there 
are no other eigenvalues with compactly supported eigenfunctions. 

Assume again that the probability measure for the i.i.d. distributed arc lengths ujj has a 
nontrivial pure point part given by (10) and consider the set 

D^i'B) := {E eR+\E = {nk/sif, Si ^ 0, Bs^/tt^ g Nq, Pi ^ 0, ke N}. 



Similarly to the analysis of Section 3 one can show that N^°°p{E) is discontinuous on the set 
C Z?^ and nowhere else. Note that for B 7^ this set is in general strictly smaller 
than defined by (1 1). Since all discontinuities of the integrated density of states are con- 
tained in N^°°P{E) this imphes that generically (if !Bsj/7r is not integer) the discontinuities 
disappear under the perturbation by a magnetic field. 

5. The Integrated Density of States and Scattering Amplitudes 

To proceed further with the analysis of the Laplacian — A(a;) we will use some results 
from scattering theory. Let Tm,n with m, n G Z, —m < n be the graph consisting of 
n + m + 1 loops Lj of arc length 00 j, j = —m, . . . , n joined symmetrically by the intervals 
of unit length and of two semi-lines attached to the loops L_„i and L„. With this graph we 
associate the Hilbert space ■Km,n = where JC*^^* = ^^(0, oo)®L^{0, 00) and 



LHLo), 



[L\lj)®L\Lj)] , m,n^O, 

=-m+l 



m 



n 



0. 



By — A^^„(a;) we denote minus the Laplacian acting on 'Km,n with the boundary con- 
ditions (4). We consider the scattering matrix and the spectral shift function for the pair of 
operators (— Aj„ „(c<j), —A) where A is a usual Laplacian on L^(]R). Although these opera- 
tors act in different Hilbert spaces, the scattering matrix as well as the spectral shift function 
can be constructed in this case (see, e.g., [23] and the A). 

Identifying in a natural way L^(M) and "K^^^ we define the isometric identification opera- 
tor a : L2(R) ^ IK such that Ran^ = Obviously, I - 3*3 = and I - 33* = ^^oxt 
with P^ext being the projection in 'Km,n onto the subspace J£'^^*. It is easy to check that 
the conditions (30) in the A are fulfilled for, e.g., k = 1. Thus the spectral shift function 
^rn,n{E', ^) '■= S,{E; —Am,n{^), "A; 3) cxists and satisfies the trace formula (31). The con- 
dition ^jn,n{E;uj) = for < fixes the spectral shift function uniquely. The scattering 
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matrix 



is defined as in Section 2 of [16]. 

The operator — Ao^ol"^) on the graph ro,o consisting of a single loop and two half-lines 
was considered in Example 3.2 in [14], where it was shown that the transmission and reflec- 
tion amplitudes are given by 

(15) To,o{E;io) = ——y= , Ro,o{E-u) = Lo,oiE;io) = - ^. 

The operator — Ao,o('^) has infinitely many eigenvalues {n'^n'^/ojQ, n G N} imbedded in 
the absolutely continuous spectrum. At those energies the reflection coefficients vanish and 
ro,o(7r2n2/a2;L.) = (-l)«+i. 

The spectral shift function ^o,o(^; can be calculated expUcitly. Indeed, by the chain 
rule for the spectral shift function we have 

CoAE; ^) = C{E; -Ao,o(u;), (-A)©(-A'°°p)) + ^{E; (-A)©(-A'°°p), -A; 3), 

where A'°°p denotes the (self-adjoint) Laplace operator on JCjj^^ = L'^{[0, loo])®L'^{[0, loq]) 
with the boundary conditions 

V;(+)(0) = ^(-)(0), ij^+y{0) + ij^-y{0) = 0, 

^^+\uo) = ^^-\uo), ^(+)'(a;o) + V^"^'(a^o) = 0. 

Combining the trace formula (31) with Birman-Krein Theorem (32) (both with with d = I) 
we obtain 

e(£;;-Ao,o(a;),(-A)©(-A'-P)) = -l0o,o(£^;u;) 

with 

ct>o,oiE-u;) := ^logdetS(£;;-Ao,o(a;),(-A)©(-Ai°°P);/) 
= Arctan tan(-\/^a;o)^ , 

where Arctan is chosen such that x h- Arctan(tan x) is continuous for all x G M and 

Arctan(O) = 0. In particular, (/>o,o(7r^P/a^; a;) = irk and (?!)o,o(7r^(A; -|- l/2)^/a^;a;) = 
Tr{k + 1/2) for all /c G No and all G il. Furthermore, 

^{E; (-A)©(-A^°°P), -A; 3) = ^{E; (-A)©(-A1°°p), -A©0; /) 

and, thus, by Lemma 3.1 in [13] equals minus the eigenvalue counting function for the 
operator -A^°°p, 

^E; (- A)©(- A^-P), - A; d) = - \VEuoM . 

Therefore, we obtain 

^0 (-2';'^) = ——Arctan ( ^ tan(v^a;o) ] — [a/Ewq/tt]. 
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In a similar way one can calculate the spectral shift function for arbitrary 

integers m and n such that —m<n, 

1 " 

j=-m 

where (f)jn^n{E; u>) is the scattering phase for the pair of operators (— Aj„ „(a;), — A). 

Denoting by No{E) = ^/E/tt the integrated density of states for the Laplacian — A on 
L^(M), by results of [13] (Theorem 4.1, equation (4.4), and Theorem 4.4) we obtain that 

(16) N{E)=No{E)- lim ^"'"^^'^^ 

m,n^oo n + m + 1 

almost surely. Although the results of [13] are formulated and proven for Schrodinger oper- 
ators on the line, all proofs extend verbatim to the present context. 



stetig_N.png 



Let 



Fig. 2. The integrated density of states N{E) for the case of uniformly 
continuously distributed ojo £ [1/2,3/2]. The horizontal plateau corre- 
sponds to a spectral gap. 



To,o{E;u>) To,o(E;w)* 
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and 



n . /o 



e+ = ( ) ~ V 1 



By the arguments presented in [13] from equation (16) we obtain the following theorem. 
Theorem 5.1. For E > the integrated density of states N{E) is given by 

(18) N(E) = T- lim '^^'''^^''^^ + iv'°°P(£;), 

TT m,n-»oo m + n + 1 

where 

(19) vi^{E;cj):=arg/e±, f[ A^.{E)e±\ , 

\ j=—m I 

(•, •) denotes the inner product in and N^°°P is given by (8). The choice of the argument 
is uniquely fixed by the condition 

for arbitrary k E 7. satisfying —m<k and —k<n. 

In equation (19) and below the product Yl is to be understood in the ordered sense. Theo- 
rem 5.1 provides an algorithm for numerical calculations of the integrated density of states. 
Figures 2 and 3 show examples of such calculations. Some other numerical results for 
Schrodinger operators on the line are presented in [12]. 

By means of Theorem 5.1 one can obtain (see [15]) a two-sided estimate on the integrated 
density of states N{E) 



(20) 



~ 1 /"°° /5 \ 

N{E) - No{E) J Arctan i- tan{u;oVE)\ d>f(a;o) 



1 

< -. 
- 2 



Comparing (18) with equations (6) and (7) we obtain that 

N{E) = ^^ lim ^-'^^^^ 



TT m,n-^oo m + n + 1 

for almost all oj e fl. A simple calculation now leads to the following representation for 

N{Ey. 

(21) N{E)=No{E)-- lim fV^ii:^, 

where (t)m,n{E; 00) is the scattering phase for the pair of operators {—^ra,n{E\ u), —A). In 
the next section we use equation (21) to prove the Thouless formula in the present context. 

Now we are in position to prove Proposition 2.1. Let Ajn,n{E;u;) denote the transfer 
matrix, 



(22) Am,n(^;w) 



1 Rm,n{E;ui) ' 

JE^ T„,„(£;;a;) 

.Tm,n{E;ui) Tm,„{E;u})* , 
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Bernoulli JSI.png 



Fig. 3. The integrated density of states N{E) for the Bernoulli distribu- 
tion X = i(52 + ^Se, i.e., ojj G {2,6} with equal probability. Horizontal 
plateaus correspond to spectral gaps, vertical strokes represent discontinu- 
ities. 

Proof of Proposition 2.1. Note that for any A; G Z is necessarily of the form a^e^'^^ + 

hke~^^^- Therefore, if ^ I /^^ ^ on an open subset of some interval Ij , then |aj|-|-|6j| ^ 0. 
For k ^ j the coefficients a^, are determined by the equation (see, e.g., [16]) 

(23) (;')=A_»,(i;;.)(«j), 

where det A^^jiE, uj) = 1. Assume that = 6fc = for some A; G Z. Then from (23) it 
follows that Uj = bj = 0, which is a contradiction. □ 

6. The Lyapunov Exponent 

Following convention we define the Lyapunov exponent for the Random Necklace Model 
as the exponential growth rate of the norm of the transfer matrix, 

(24) jiE)= lim ^o^W^rnAE^l 

m,n->oo n + m + 1 
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It is a general fact that for every E > this limit exists almost surely and is nonnegative. 
Moreover, (see Theorems 5.1 and 5.3 in [13], also cf. [19], [20], [21]) 



(25) 



7(E) = - lim 

m,n—>-oo 



log\\Tm,n{E;uj) 
n + m+ 1 



where Tm,n{E; to) is the transmission coefficient corresponding to the Laplacian — „(a;). 
Also, (24) can be rewritten as 



(26) 



7(E) 



lim 



m,n-*oo m + n + 1 



= lim 



m,n-^oo m + n+1 



log 



log 



n A-.(^) 



where the matrices A^^^. (E) are defined by (17). 

Kotani's Theorem, which states that the Lyapunov exponent of every Schrodinger opera- 
tor on the line with non-deterministic potential is almost everywhere positive (see, e.g., [3]), 
does not apply directly to the model considered here. Therefore, to prove that 7(-B) is pos- 
itive for almost all E > in the case when supp k contains at least one non-isolated point 
we refer to Theorem 5.6 of [13]. By this result for E > the Lyapunov exponent vanishes 



on the set {E 



n e N} and nowhere else. Also, for E = the matrix A^„ (E) 



equals the unit matrix and hence 7(0) = 0. As an illustration we have computed the the 
Lyapunov exponent for uniformly continuously distributed u>/s on the interval [1/2, 3/2], 
i.e., for dK{E) = X[i/2,3/2](^) (see Figure 4). 

The case when the measure x is purely discrete (and thus is a finite convex combination 
of Dirac measures) is not covered by Theorem 5.6 in [13]. First, consider the Bernoulli 
distribution x = pSg^ -|- (1 — p)6s2, < p < 1 since in this case we may invoke the 
recent results of Damanik, Sims, and Stolz [5] (see also [6]). To apply this result we need 
to introduce the scattering amplitudes T^^i'^"), R^^^'^°\ for — Ao,o(si) relative to 

the "background" operator — Ao,o('So)- (As discussed above the fact that these operators 
act in different Hilbert spaces is not relevant). A simple calculation shows that they can be 
determined from the relation 




_R(^l.ao)(E) ' 

"t(''1'*o)(_E) 



Ao,o(E;si) Ao,o(E;so) 



where Ao,o is defined in (22). In particular, we obtain 



j2iso-/B^ 



such that the reflection coefficients vanish if and only if y/E(si — sq) /tt is an integer. There- 
fore, by applying Theorem 1 of [5] we conclude that the Lyapunov exponent vanishes on 
the set 



(27) S{so, si):={E=l — ] ,keN}u{E 



irk 
2 



nk 



Sl - So 



ken 



and nowhere else. 



A RANDOM NECKLACE MODEL 



13 



stetig_gam.png 



Fig. 4. The Lyapunov exponent for the case of uniformly continuously 
distributed a;o e [1/2,3/2]. 

Let now x be an arbitrary discrete measure given by (10) with a finite number of non- 
trivial terms. As noted in [5] the set of zeroes of the Lyapunov exponent is contained in the 
union 

U S{s,s') 

s^s' 
Sjs'esupp X 

of the sets (27). The Lyapunov exponent is strictly positive for all E away from this discrete 
set. 

Results of numerical computations of the Lyapunov exponent using (26) for the BemoulU 
distribution with x = ^62 + ^Sq, i.e., ujj G {2, 6} with equal probability, are presented in 
Figure 5. We mention two properties of 7(£'): 

1. The Lyapunov exponent is periodic in \/^ with period vr, i.e. 

7((\/^ + 7rA;)2) =7(E), ^ > 0, fc G N. 
This follows immediately from the fact that for even s 

AsiiVE + ^f) = As{E). 

2. The Lyapunov exponent is reflection symmetric, 

(28) 7((A;7r - Ve^) = -^{E), E > G N, k> /E/tt, 
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Fig. 5. The Lyapunov exponent for the BernouUi distribution x = + 
^(^6, i.e., (jjj G {2, 6} with equal probability. 

i.e., the points /c G N are the axes of reflection symmetry (on the scale of Ve/tt). To 
prove this we note that 

Since 

up to a sign the product nj=-n ^^^j ((^^ ~ V^)^) equals 

(rr)nM-)(rr)- 

\ / j=—n ^ ' 

Thus, equality (28) follows from (26). 

Using equation (26) one can also analyze the periodic case Uj = loq for all j G Z. In this 
case the spectrum of — A(u;) consists of the absolutely continuous part and the eigenvalues 
Ek = TT^k^/uQ, k G N of infinite multiplicity. The absolutely continuous spectrum has a 
band structure such that E G speCa(,(— A(c<;)) if and only if Hill's discriminant (cf. equation 
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(8) in [2]) 
H{E) 



1 



cos{^/E{uo + 1)) - - cos{Ve{u}o - 1)) 



\To,o{E;u)\ 4"™^' ^"^ ' 4 

satisfies the inequality \H{E)\ < 2. From this and the fact that |ro,o(7r^n^/a^;a;)| = 
1 it follows immediately that all eigenvalues are imbedded in the absolutely continuous 
spectram or he at the edges of the spectral bands. 

7. Continuity of N{E) 

For almost all > the Lyapunov exponent satisfies the Thouless formula 

\-E 



(29) 



7(E) = a + / log 



A 



dN{\), 



where a is some positive number. The existence of the integral on the r.h.s. is guaranteed by 
the estimate (20) and Lemma 1 1.7 in [22]. We emphasize that N^°^p{E) does not enter this 
formula. 

Before we proceed, we discuss the implications of (29) to the continuity properties of 
N{E). By a modification of an argument due to Craig and Simon [4] (see Theorem 11.9 in 
[22]) the positivity of the Lyapunov exponent imphes the log-Holder continuity of N{E), 
that is, the inequality 

NiE2)-N{Ei) <C\log\E2-Ei\\-'^ 

for arbitrary sufficiently small intervals [Ei,E2]. Moreover, using the arguments (in a 
shghtly modified form) of Damanik, Sims, and Stolz from [6] one proves that N{E) is 
actually Holder continuous, i.e., there is a number < fi < 1 such that 

N{E2) - N{Ei) < C\E2 - Eil^" 

for arbitrary sufficiently small intervals [Ei, E2] not containing the points where the Lya- 
punov exponent vanishes. Both these properties hold for general distribution measures x. 

We only sketch the proof of (29), but the interested reader can easily fill in the details. 
We closely follow the line of arguments given in [13]. 

Let f{z) be an analytic function in open cut plane Co = C \ [0, 00), continuous in the 
closure of Co and satisfying \f{z)\ < c\f\z\ for all z G Co- Moreover, we assume that 
+ ^0) = f{E — iO) is continuously differentiable for all E > 0. Using the Cauchy 
integral formula it is easy to show that 



Refiz) 



Re/(i) + - 

TT 



log 



A 



A 



dIm/(A + iO). 



From equations (15) and by the factorization rule for the matrices (22) (see [16]) it follows 
that one can take f{z) = log Tm,n{z\ oj) for arbitrary m, n, and 00 e O,, thus, obtaining 



1 



log \Tm,niz;Uj)\ = log \Tm,n{i; + - log 



TT 



A 



A 



<i^m,n(A;a;). 



Using Lemma 1 1 .7 in [22] we conclude that for almost all > 0, all to E il, and arbitrary 
integers m, n 



log\Tm,niE;Uj)\ = log|r^,n(i;w)| + 



1 



TT 



log 



X-E 
X-i 
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Now divide both sides of this equation by n + m + 1 and consider the limit m, n oo. 
By (25) the l.h.s. converges almost surely to —j{E). The signed measures (n + m + 

l)~^d (pni.n(E\uj) converge vaguely to Tr{dNo{E) — dN{E)). Thus, again by Theorem 
11.7 in [22] there are subsequences m^, such that 



lim — 

fc^oo vr 



log 



A 



ti'/'mfc,«fc(A;u;) 
mk + Uk + l 



/ log 


A-^ 




X-i 



for almost all E. Noting that 



where jo{z) 
obtain (29). 





A-^ 


/ log 


/k 


X-i 



{dNo{X) - dN{X)) 



V2 



dNo{X)=jo{E)-jo(^) = -^, 
I Re \/^| is the Lyapunov exponent of the Laplacian — A on 



we 



Appendix A. The Spectral Shift Function 

Here we briefly collect some facts from the theory of the spectral shift function in the 
case where the operators involved act in different Hilbert spaces. For a comprehensive 
presentation we refer the reader to the book [23]. Consider two (possibly unbounded) self- 
adjoint operators Tq > I and T > I acting in Hilbert spaces "Kq and J{, respectively, and a 
bounded operator 3 '■ "Kq ^ 'K. Suppose that the operators 

(30) {d*d-I)TQ'^ \ are trace class. 

T-\33*-I) J 

Under these conditions there exists a spectral shift function ^ {E; T, Tq ; 3) for which the 
trace formula 

(31) tr [ct>{T) - a,^(ro)a*] +tr [{3*3 - mn)] = [ ^{E-T,To;3)<i>'{E)dE 

JR 

holds, where ^ is an arbitrary bounded continuously differentiable complex- valued function. 
The relation to the scattering matrix is given by the Birman-Krein theorem 

(32) det S{E; T, Tq; 3) = exp{-27r^e(£^; T, Tq; a)}, 
where S{E; T, To;3) is the scattering matrix for the operators (T, Tq). 
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